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Pj'opagation through a gyromagnetic medium in a direction normal to a 
uniform magnetizing field is considered. Geometrical arrangements which 
make this propagation non-reciprocal are described. A few illustrative 
examples are discussed briefly. The non-reciprocal helix, of importance in 
traveling wave tube work, is treated at length. 

1. INTRODUCTION 

1.1. General Remarks about Non-Reciprocal Propagation 

Part I of this paper began with a brief discussion of some of the micro- 
wave properties of two gyromagnetic media; the gas discharge plasma 
and the ferrite. The remainder of Part I was devoted to the analysis of 
the mode spectrum in a cylindrical waveguide filled with one of the 
media and placed in an axial magnetic field. It was demonstrated that 
the natural modes in such a guide are right- and left-circularly polarized 
waves which travel with different phase velocities. Accordingly a plane 
polarized mode, which to some approximation can be regarded as the 
sum of right and left circular modes, will, in traversing a section of the 
guide, undergo Faraday rotation, just like a plane wave in the un- 
bounded medium. It is true that the presence of the guide wall has a 
drastic effect on the course of the rotation with magnetizing field, chang- 
ing it, sometimes beyond recognition, from that prevailing in the un- 
bounded medium. Nevertheless the principle remains the same; confine- 
ment of the wave to a guide merely modifies quantitatively the Faraday 
effect for plane waves. In optics, where practically plane waves are 
almost always employed in this connection, the non-reciprocal nature 
of this effect is so familiar that it hardly requires restatement here. 
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By contrast, Part II of this paper deals with devices whose non- 
reciprocal operation depends in principle as well as in numerical detail 
on the disposition of the boundary, or, more generally, on geometrical 
configuration. These devices employ magnetizing fields transverse to 
the propagation direction. Some electromagnetic field configurations are 
unaffected by such a dc field, but, whenever the rf magnetic field in the 
case of a ferrite (or the rf electric field for a plasma) has a component 
normal to the dc magnetic field, this is no longer the case. For, now, the 
magnetization (or the charges) will be caused to precess about the dc 
field, giving extra terms in Maxwell's equations and a resultant change 
in the propagation. This change may be simply an alteration in phase 
velocity, the propagation remaining reciprocal. This is the case, for 
example, for the propagation of plane waves in an infinitely extended 
medium [Cotton-Mouton effect]. Here, since every direction of propaga- 
tion normal to the dc field is physically equivalent to any other and, in 
particular, to the opposite direction, no non-reciprocity can arise. 

For reciprocity to be preserved in the presence of the dc magnetic 
field is, however, exceptional and requires a certain amount of geo- 
metrical symmetry in the system. That non-reciprocity may be expected 
in asymmetrical systems may be foreseen if we consider a system, typical 
of those to be treated in this paper, in which all the rf fields are hide-, 
pendent of the coordinate along which the dc magnetic field is pointing. 
The relevant conducting boundaries and any interfaces between ferrite 
(plasma) and air are all surfaces parallel to the direction of propagation 
and lying in the dc magnetic field direction. Suppose the system to be 
divided into two parts by another surface of a similar kind and examine 
the surface impedance of one of the parts (which should contain some 
gyromagnetic material). If the propagation direction be reversed it is 
necessary to reverse the magnetic field to retrieve a situation in the 
part considered geometrically equivalent to the original. But, since the 
precession of the magnetization (or charges) about the dc magnetic 
field has a definite sense, the magnetic or electric current associated 
with this precession will be reversed when the dc field is reversed. Thus, 
the properties of the medium are altered and the surface impedance will 
be different for the two directions of propagation. In general, the surface 
impedance of the other part of the system will not compensate for this 
distinction between the two directions and we shall find different propa- 
gation constants for opposing directions. An exception will occur if the 
system contains a surface about which it has geometrical symmetry, for, 
then, compensation clearly takes place about this surface and the system 
is reciprocal. 
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An example of a simple non -reciprocal system is indicated in Fig. 1(a). 
Here a slab of ferrite is inserted into a rectangular waveguide parallel 
to the narrow walls and closer to one of them. Several workers have 
demonstrated that this arrangement and a similar arrangement in a 
circular waveguide are non-reciprocal for what is essentially the dom- 
inant mode. 1, 3 * When the slab is centered in the guide we have a plane 
of symmetry and the non-reciprocity vanishes. 

Another configuration of the transverse field type is represented by 
the system shown in Fig. 1(b). Here a hollow ferrite cylinder is mag- 
netized circumferentially and propagates a TE 0n -mode. It is clear that 
any arrangement of this sort, which might, in principle, include conduct- 
ing sheaths, internally or externally, or might have the ferrite extending 
to an indefinitely large or small radius, cannot have any symmetry 



WAVEGUIDE - 



WINDINGS FOR 

MAGNETIZING- 

CURRENT 




PROPAGATION 
DIRECTION 



CIRCUMFERENTIAL,-' 
MAGNETIC FIELD 

IbJ 

Fig. 1 — (a) Rectangular waveguide and ferrite 
magnetized ferrite cylinder. 

* It is expected that an article on this subject by S. E. Miller, A. G. Fox and 
M. T. Weiss will appear in a forthcoming issue of the Journal. 
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about a cylinder coaxial with the ferrite. Thus, the internal and external 
impedances of such a system at any coaxial cylinder can only com- 
pensate accidentally (perhaps at a single frequency) and non-reciprocity 

is the rule. 

It is frequently asserted without qualification that for non-reciprocity 
a further condition upon the relevant rf field is that its projection upon a 
plane normal to the magnetizing field be elliptically or circularly polar- 
ized in the limit of vanishing magnetization. The argument is based on 
the consideration, in itself correct, that the effective material constants 
are different for right- and left-circularly polarized field vectors. Suppose 
that the magnetization direction is y. Then the tensor relating B to H 
is (see Part I, Section 2.1): 







Mo 



-J« 



JK U / 

For right- and left-circular fields withH z = ±jH x , therefore, the medium 
is isotropic in the plane transverse to the dc field ^vith permeabilities 
^ _j_ K} n — K respectively. Since opposite circular polarizations accom- 
pany opposite propagation directions, (see for example, Fig. 2) the per- 
meabilities, and hence the propagation constants, are different for oppo- 
site propagation directions. It is then argued that the field must already 
be circularly, or at least elliptically polarized to start with, if non- 
reciprocal effects are to result from application of the magnetization. 
However, the argument is true only for effects of first order in the mag- 
netization. For general values of magnetization, the rf field, even if 
linearly polarized to begin with, will become elliptically polarized, and 
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Fig. 2 — Magnetic lines of force parallel to the broad side of a rectangular 
waveguide. 
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non-reciprocity will occur. It is understandable that this double function 
of the magnetization (conversion to elliptic polarization plus creation 
of differences in permeability) leads to higher order effects. However, 
inasmuch as for the ferrite and plasma the magnetization can produce 
large changes, the requirement of elliptic polarization in zero magnetic 
field cannot be regarded as essential in practice. These considerations are 
demonstrated by a simple example in Section 2.2. 

The devices considered in this paper actually are such that the electric 
or the magnetic field vector in the plane normal to the field is elliptically 
polarized even in the absence of the gyromagnetic medium. For example, 
the magnetic lines of force of the TEi mode in a rectangular waveguide 
form two sets of closed loops in a plane parallel to the wide sides (see 
Fig. 2) and repeating every wavelength. This pattern moves bodily 
down the guide with the phase velocity of the mode, so that an observer 
stationary at any point not at the center or at the narrow walls of the 
guide sees a magnetic field rotating at the signal frequency and tracing 
out a generally elliptic path. The sense of the rotation is opposite on 
opposite sides of the center plane, and depends on the propagation direc- 
tion. The conditions outlined in the previous paragraph are therefore 
satisfied; introduction of a ferrite slab magnetized as in Fig. 1(a) will 
yield first order non-reciprocal effects. 

The problems considered here are such that the electromagnetic 
fields do not vary in the direction of magnetization. Under these condi- 
tions the field can be split into a TE and TM field satisfying different 
wave equations. In general, the two fields are coupled through the bound- 
ary conditions. Most of the paper is devoted to the analysis of the non- 
reciprocal helix, a problem that has recently gained importance in 
connection with high power traveling-wave-amplifiers. 4 The conventional 
amplifier suffers from a limitation on its maximum useful gain; waves 
reflected from the output end will make the tube "sing" above a certain 
critical gain. Ferrites offer the possibility of preferentially suppressing 
these backward waves and so of increasing the permissible gain by a 
large amount.* In section 2.3 the "flat" helix (one of infinite radius) is 
considered. For the slow waves employed in practice a rather complete 
treatment is possible in this case of planar geometry. In Section 3 the 
cylindrical helix is treated. Inasmuch as the solutions involve functions 
for which no extensive tables exist, the treatment has to be more sketchy. 

* More specifically, in high-power traveling-wave tubes the large beam cur- 
rent employed may be above the critical value required fo r backward wave os- 
cillations due to spatial harmonics of the helix structure. In such cases the larger 
attenuation of backward waves will permit a higher beam current and therefore 
stable amplification to higher power levels. 
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Thus the loss is neglected and only the non-reciprocal phase charac- 
teristic is considered. The losses have then to be determined approxi- 
mately by differentiation of the phase characteristics. 

A few further problems were considered as illustrations of the general 
principles. One case, that of the plasma filling the space above an im- 
pedance sheet can actually be solved analytically and provides a par- 
ticularly clear demonstration of the non-reciprocity. The case of the 
rectangular waveguide with a ferrite slab has already been considered 
extensively elsewhere, and only the results for a thin slab are given here. 
A problem with cylindrical symmetry is taken up in Section 3.3: a 
cylindrical waveguide fitted with a circumferentially magnetized cylinder 
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Fig. 3 — p a versus a for various p. 
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of ferrite. Again the discussion had to be sketchy in view of the scarcity 
of information on the functions that solve the problem. 

2. PLANAR GEOMETRY 

2.1. Fields and Impedances 

In this section we consider planar transverse field problems which are 
characterized by the following conditions. The dc magnetic field is of 
uniform strength H within the gyromagnetic medium and points along 
the ?/-axis. All rf field components are independent of the y coordinate. 
We discuss the ferrite case first, then indicate how the results are to be 
translated for the plasma. 

For the orientation of the dc magnetic field which has been chosen the 
permeability matrix is of the form: 



3* 







Mo 





->- 



M and k are, in general, even and odd functions of H ; the permeability 
of unmagnetized ferrite is taken to be no as in free space. Following the 
procedure of Part I we shall assume specifically for n and k the formulae 
given by Polder's treatment of the dynamics of the medium. Thus, we 
have the expressions (for the case of no loss): 



H 1 — pa — a" 
Mo 



1 - (T 2 



K p 

Mo 
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K p 

Ph — = ; , 

/x 1 — pa — a 1 



where a is the ratio of the precession frequency, L— # , to the signal 



7_ 
2tt 



frequency and p is the ratio of a frequency, L-— ■ M /no , associated with 

2r 

the saturation magnetization, M , to the signal frequency. It should be 

noted that p and a have always the same sign. The behavior of n and k 

as functions of a was shown in Fig. 1(a) and (b) of Part I. p H is shown as 

a function of a in Fig. 3. The dielectric constant of the ferrite is taken to 

be €. For reasons given in Part I, | p \ is assumed less than unity. 
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When the condition, d/(dy) = 0, is put into Maxwell's equations the 
latter are found to be separable into two sets: 

_ ?*& - 4**., (2a) 

dz 

S-jWfc, (2b) 

ox 

W. _ 9B. _ _ ^ oH|/ (2o) 

and 

_ *S _ - jwL^, - /*ETJ , (3a) 

oz 

a ^= -jMfjkttp+MffJ, (3b) 

d A*- d JL> = jcoeE v . (3o) 

It is to be stressed that such a separability is possible only when the rf 
fields do not vary along the dc magnetic field. The sets of equations (2) 
and (3) correspond to the separate equations for H z and E z which arise 
from (13) of Part I when is there set equal to zero. The first set de- 
scribes a TM field of the familiar type, whose propagation through the 
medium is unaffected by the magnetic field. The second set describes 
a TE field whose components, because of the presence of k, are con- 
nected by different relations from those which exist in an unmagnetized 
medium. The separability of the two fields is equivalent to saying that 
they are not coupled by the medium itself, but they may, of course, be 
coupled at the boundaries. 
We may write (3a) and (3b) in the form 

- H(X 2 -^)H x = j^-^, (4a) 

- M S-* 2 )H,= ^ + j* d -§?, (4b) 



and upon eliminating H x and H z , the wave equation for E y is found to be : 

**+**+*£^*M.-0, (5) 

dx 2 dz 2 n 

where E v (and also the H's) are evidently propagated in the ferrite as 
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though the latter had an effective permeability, (p — k)/h = /i eff . This 
permeability may assume any value from + <x to — <*> as may be seen 
by writing it in terms of a and p. From the Polder formulae, in fact, 



/ 2 2 W 1 - (p + a-) 2 

Meff - (.M — K ;/M = MO 



1 — pa — a 2 



(6) 



As it should be, this is an even function of magnetic field. Meff/Mo decreases 
from 1 — p to as | a \ rises from to 1 — | p \ ; it decreases from to 

— » as | a | runs from 1 — | p \ to \/l + p 2 /4 — I V I/ 2 a nd finally 
decreases from » to 1 as | <r | increases indefinitely above \/l + p 2 /4 

— | p |/2. Its behavior is indicated in Fig. 4. It should be recalled from 
Part I that "a = 0" is an abbreviation for the very small magnetic field 
necessary to saturate the ferrite. 

A brief examination of the propagation of plane waves shows the more 
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significant features of transmission in this medium. Since no direction 
in the x-z plane can be a preferred one, the plane wave may be assumed 
to travel in the z-direction with variation, e"*\ Then, from equations 
(4) and (5) 



and 



J = co en eS 



H* T~2 2\ ^11 > 

W(/» a — K z ) 

C0(]U J — K £ ) 



Since n eS is negative between \a\ = 1 — \p\ and | a \ = \/l + p 2 /4 
— | p |/2, the medium is cut off for plane waves in this range of magnetic 
field (at a fixed signal frequency). H is elliptically polarized when the 
medium is magnetized and \H,\/\H X \ = \ k |/| m I = I Pm |- We may 
also put 

H x + jH. = 



H x - jH z = 



w(/x — k) ' 
co(/x + k) ' 



But | H x + jH z | and | H x — jH t \ are proportional to the amplitudes of 
the left-handed and right-handed circularly polarized components of 
the magnetic field. The medium may thus be considered to exhibit the 
permeability, 



/J. — K = HO I i — 



1 - 

for left-handed components and the permeability, 
M+K = M0 (l+_P_j 

for right-handed components. The effective permeability is essentially 
a parallel combination of these two permeabilities and the medium may 
propagate (n eB > 0) even when m — k is negative, as will happen for 

Vl + P 2 /4 - I V 1/2 < <r < 1. 

Since the medium itself has no non-reciprocal properties it is clear 
that if the latter are to arise they must do so as a result of interaction 
between the medium and its surroundings. The boundaries of the ferrite 
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at which matching 'will be necessary are surfaces parallel to the y-axis 
and here we need an expression for H tan JE y where H taaK is a tangential 
magnetic field at the surface. For the moment we will consider the 
admittance looking into the ferrite and take tangential components in 
the counterclockwise sense. From equation (4), then, 

. \_BEy 1 8E y 

H ta ng _ E y dv E y da (7) 

~J^' " 0)(M 2 - K 2 ) W( M 2 - K 2 ) ' 

where d/(dv) is a normal derivative (outward) and d/(da), a tangential 
derivative. It is possible, although by no means essential, to interpret 
the terms of (7) in the following fashion: the first term is just the admit- 
tance of a normal TE mode propagating in the interior of the ferrite 
(which is to have the permeability, n eB ) ; the second term is to be ascribed 
to an independent surface current, 

dE v 

— K 

da 



0>(/i 2 ~ K 2 ) ' 

Using this picture one may see how non-reciprocity arises in a simple 
case. If the ferrite be bounded by the planes x = Xi and.r = x 2 , and the 
z-variation is of the form e - **, the admittances due to the surface cur- 
rents are +jk@/u(h 2 — k). If /3 reverses its sign the surface currents are 
interchanged. If now the external admittances on the two sides of the 
slab are unequal (and, of course, themselves reciprocal) for given values 
of w and \&\, there is no reason why /3 and —0 should simultaneously 
solve the matching problem. 

Almost all the above considerations may be taken over to the case of 
the plasma. Here the TE fields will be undisturbed by the magnetic field 
(but the dielectric constant is altered by the presence of the charge 
from its free space value) . Equations (4) and (5) are now replaced by 

M S- v *) El = -.**£ + **%!, (8a) 

M^-^^T'+e^, (8b) 

az ox 

■^r + ^r +WM0 ^— H "-° (9) 

for the TM fields. Here t and -q are the diagonal and off -diagonal terms 
of the dielectric matrix which is of the same form 
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as the permeability matrix of the ferrite. The equations of motion for 
the plasma* lead to the expressions 

ei = € (1 — q ), 



e = e 1 + 



- 1 



V = Co 



2 

qa 



(10) 



<r 2 - 1' 

2 



_ €" - if _ (1 - g ) - 0- 

€eff e €0 (1 - q*) - o 2 ' 

where a is now the ratio of the cyclotron resonance frequency; 

1 I e I mo 



2tt m 



H, 



in a field H , to the applied frequency and q is the ratio of the plasma 
frequency to applied frequency; e e ff behaves with magnetic field in much 
the same way as n eS . It is a constantly decrea sing fun ction of a and is 
negative between | a- 1 = 1 - q and | a \ = Vl - q 2 , g° in g to infinity 
at the latter value. The left and right handed dielectric constants, 
e - v and e + t? are given by e [l - q/(l - o-)]and« [l - q 2 /(l + <r)]. 

2.2. Examples of Non-reciprocal Systems 

We now discuss briefly three examples of non-reciprocal systems as 
illustrations of particular points. As an example of a system which can 
be analyzed very easily and completely we consider a plasma occupying 
the region, x > and bounded at x = by a sheet of constant impedance. 
This impedance is to depend upo n freq uency but not on the propagation 
constant. It will be written as jVmoAoZO) where /i and e are free space 
values. A practical realization of such a sheet might consist of a very 
large number of similar fins of negligible thickness and separation, 
parallel to the ?/-axis and attached normally to a conducting plane x = 
constant. The fields between separate fins are uncoupled and E , is u ni- 
form between fins. For such an arrangement, Z(o>) = tan o>VeoM'o.ro » 

* See Section 2.2 of Part I. 
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where .r is the depth of the fins. If the z-variation of the fields is e~' p ' 
and the waves are guided, the x-dependence in the plasma must be 
as exp — \//3 2 — or/zoceff x - From equation (8b) 

j<o(e — r] )E, = [v6 - eV/3 2 - ovWeff ]H V . 

Matching at x = gives 

?7/3 - e V/3 2 - o) 2 /xoe e ff _ . Ao 7 / ,\ 
Me 2 - rf) V € 

This yields 

\fi - ijco i/S Z(«)T = a>W + "- e 2 Z 2 (co) 
r co J fo 

or 

P 



Moto 



= I Z(«) ± 4/1 + - 9 Z 2 (o>) . (11) 

The non-reciprocity is clearly exhibited, since the two values of /3 are 
not equal and opposite. The solution (11) is valid only if /3 2 > co"W e ff , 
corresponding to guided waves. 

In the second example we assume that the region between conducting 
planes at x = and x = Xo is filled by a plasma. When no magnetic field 
is present E z is supposed to vanish and E x is uniform across the gap. 
The unperturbed field is then plane polarized (TEM). The magnetic 
field is now applied parallel to the y axis to that part of the gap between 
x = and x = Xi . E t in the magnetized plasma is now given by 



E, = E sin Vco' 2 Moeeff — /3 2 z 

since it must vanish at x = 0. The z variation is again exp — jflz. H y 
may be found from equation (8b) and is 



H„ = 



3<a ° M [nfi sin V«Vn€eir - /3 2 x 



— e V^Weff — /3 2 COS VwVoCeff — & 1 A- 

The admittance of the magnetized section at x = X\ , is thus, 



r , b/3 - € Vw'-Weff - /3 2 COt VwVofefl = /3 2 Zl], 

and, analogously, that of the unmagnetized part is 

Jf = ..?... JU o i [ _ei Vw 2 Mo€i - 2 COt VoWi - /3 2 (Xo - Ii)], 



, 
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where ei = e (l — q). Suppose now that the applied field is weak, so that 
« 2 juo<: ~ co 2 /k>€i ~ /3 2 - Then, the cotangents may be expanded and by 
equating admittances one obtains 

jg f" a _ ±\ = Ml 

coVoe - /3 2 |_ zj («Wi - /3 2 )(z - xj ' 



or 



f =1 + 



-«i + Uj8 - jM (»o - »i) 



W"/io€i \ ( a e 



(12) 



Since e — ej is of the second order in o- this equation may be solved by 
substituting the unperturbed value of /3, ± w\/mo«] » iu the right hand 
side. It is clear that although the system is non-reciprocal, it will be 
so only to third order in a. This system, therefore, illustrates the fact 
pointed out in Section 1.1 that even when the fields are plane polarized 
in the absence of a dc magnetic field, non-reciprocity may arise, although 
it may be very small in weak fields. 

The third example to be considered is one which has been referred to 
in Section 1.1, namely that in which a strip of ferrite is placed across 
the short dimension of rectangular wave guide, see Fig. 1(a). In view 
of the fact that this problem has been discussed with great thoroughness 
by Lax, Button and Roth, we shall, after deriving the characteristic 
equation, consider only the case of a very thin strip. Let the thickness of 
the strip be 2x and the distance from its two faces to the nearest guide 
wall be Xi and Xz respectively. The admittances at the two faces are then 



and 



-2- cot V/V - /3 2 xi 

o>Mo 



—? cot V/V - /3 2 Xi 

CCfXo 



respectively, where /3 2 = w 2 juoeo • Inside the ferrite 



jg- -t«V«*-rt- -to 8 -A 



and 



OX 
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Thus, immediately within the ferrite, 

1 dE u (H,\ 

Iliy OX \iiy/external 

If the two faces of the ferrite are x = — xo and x = xo we then have 

/ i dEA = gj. cot v ^r^ Xl _ PBfi = A , 

\Ly dX /l=-X HO 

and 

(w ir) = '— cot Vpo^J* x, - prf = B. 
\hy dX /x-x Ho 

If we write 

i ^ = V ^^ 2 tan (V^7^7 2 *), 

and make use of the boundary conditions we obtain 



^^-^■ ti?-^ - 



(13) 



The non-reciprocity is clearly contained in the odd power of /3 in A 
and B. 

For small thickness we replace the tangent by its argument, and, sub- 
stituting for A and B on one side of the equation, obtain 

^5 [cot V/V - /3 2 Xi + cot V/V - /3 2 xj = 2rr [/3/ 2 - /3 2 - A£], 

Mo 

or 

sin V/V - /3 2 (»i + a;..) 

(14) 

= 2xo — sin V/3o 2 - /3 2 zi sin V/3 2 - /3 2 xa[0/ 2 - /3 2 - AB]. 

Meff 

Since the guide is almost empty, we may write 

V/So* - /8»(xj + a*) = * - 5, 
where 5 is small. Or, 

/3i(«i + 2:2) 2 
where ft 2 = /3 2 — x 2 /(xi + Z2) 2 . 
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Writing 8 = ft in the right hand side of equation (14) and noting that if 



Vft 2 - P? Xi = 



irxi 



Xi + x 2 
then Vft 2 - ft 2 x 2 = t - 6, we have 

8 = 2X0-^- sin 2 6 L 2 - ft 2 - fej! cot - P*ftY^ cot - ftft)] , 
= 2xo ^ [(ft 2 - (1 + p/)/3i 2 ) sin 2 - (fisj cos 2 (15) 

+ 2 ^5 p H /3 x cos sin . 
Mo J 

The non-reciprocal part of 8 is thus 4ikco(zi + x 2 )~ 2 p H sin 20. This has 
a maximum value for = r/4 or 3tt/4 and, hence, cd = (xi + z 2 )/4 or 
3(zi + rc 2 )/4. This result may be understood qualitatively by con- 
sidering the fields in the guide before the ferrite is inserted. We then have 
E v = E sin (vx/a) where a = x x + x 2 and consequently, 

. T 

J — 

8 . 7rx , TT a tX 

Hx = —— sin — and H s = — cos — . 

The amplitudes of the left- and right-handed components of circular 
polarization at a point are then proportional to 

. tx t tx , . tx , t tx 

— 3 S in — — - cos — and — 8 sin — + - cos — . 

a a a a a a 

The difference in the squares of these amplitudes is 26 (t/ a) sin (tx/o,) 
cos (tx/u) and this is proportional to the difference between the energy 
stored at x in the left-handed wave and in the right-handed wave. It 
is plausible that this should be a measure of the non-reciprocal effect 
produced by a thin piece of ferrite at x. 

2.3. The Plane Helix 

In dealing with transverse field problems with cylindrical geometry 
we shall consider non-reciprocal propagation along a helix which is 
surrounded by circumferentially magnetized ferrite. The analysis of this 
problem is rather cumbersome and it is advantageous to study first an 
analogous plane problem. The simplicity thus gained allows us to 
examine somewhat more complicated problems. The "plane helix," 
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to be treated here, is a sheet of negligible thickness, lying in the plane 
x = 0, which conducts only in a single direction making an angle \f/ 
with the ?/-axis. In this direction it will be supposed lossless. In addition 
we assume that the regions, x < and < x < x are empty, while the 
space xi > x is filled with ferrite. As usual the magnetic field is along 
the ?/-axis and the fields are independent of y. The problem is clearly 
the limit for very large radius of that of an empty, helically-conducting 
cylinder, surrounded by an infinitely thick shell of ferrite, circumferen- 
tially magnetized, the whole system carrying fields with no angular 
variation. 

We first consider the boundary conditions for the plane helix, after 
noting that it is evident that both TE and TM fields will be required. 
The tangential electric field on either side of the sheet must necessarily 
be at right angles to the direction of conduction since the conductivity 
is infinite. Further, the tangential electric field must be continuous through 
the sheet. Hence, if the field normal to the direction of conduction is 
E (omitting here and elsewhere the factor e"**), we have 

E a + = E z " = E cos ^, 

E y + = E v ~ = -E sin ip, 

where the symbols + and — refer to x > and x < respectively. 
Again, since current cannot flow normal to the direction of conduction, 
the tangential magnetic field along the latter must be continuous 
through the sheet or 

(H, + - H z ~) sin 4> + (H y + - H y ~) cos f = 0. 

The boundary conditions may be combined into a single equation, by 
introducing admittances, in the form 

H z H Z \ . 2 , I Hy Hy \ 2 • /, n\ 

— =- ) sin f = I ■==- — — ^ I cos \f/. (16) 



Ey+ Ey-J ' \E+ E Z -. 

The left-hand side refers to the TE fields and the right to TM fields. 
In the empty regions surrounding the sheet 



j~ 2 - (0 2 - w 2 e /xo) H v 



h - ^ - **> 



Hy = 0, 



F 1 9Hy 

JW€ OX 



and 
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dx 2 



2 )] Ey = 0, 



-1 dE y 

tiz — -. -r— , 

JCO/iQ oX 



with /3 2 = w 2 €oMo • If the w aves are to be guided, /3 2 > /3 2 , and, then, for 
x < 0, we have d/dx a V/3 2 - /3 2 . Thus 

#.- V/3 2 - /3o 2 ' 
and 



V/3 2 - /So 2 



jwyxo 



If the admittances at the surface of the f errite are Hj/Ej and H//E v f , 
then H v + /E z + and H z + /E y + are given by the impedance transformation: 

V<y . ~ A' f£ - tanh V^ft" 2 *o 
; we o JU e El M 



5l = 

E+ V/3 2 ~ /3o 2 



1 _v^^^ tanhV __ 2 



jcoe 



rc 



and 



JW/io 



# 



Hi = V/3 2 -/3 Q 2 . V/3 2 -/3o 2 .&V 
#„+ jcoeo , — J«Mo jff/ 



3T 2 JT, - tanh Vj8» - ft* *» 



1 - 



_V/3 2 - /3o 2 «/. ' 



•tanh V/3 2 - /3 2 xo 



Within the ferrite the TM-fields fa ll off as exp - V/3 2 - coVo^ and the 
TE-fields as exp - V/3 2 - uPepefP- We then have 



*/_ 



— jcoe 



&' V/3 2 - «W 



and 



ffi 



J«Mefl 



E y V/3 2 ~ «W - fa/? 

These results may now be collected and substituted in equations (16). 
The equation of condition so obtained is 

A + 1 



(/3 2 - /3o 2 ) 



A + tanh V/3 2 - /3 2 Zo 

-ft'^' i-BtJl.'v^-MV (17) 
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where 

A = Meff Vp ~ 0O 2 



Mo V/3 2 - W 2 €Meff - 0Pff 

and 

e V0 2 - 0o 2 



B - - 



«o V0 2 — w 2 



ep.n 



We shall assume that we are dealing with slow waves (0 large). This 
is the case of greatest practical interest and is usually ensured by making 
cot \J/ large. Assuming that the waves are slow we simplify the equation 
(12) and find certain values for 0. We may then ascertain in what ranges 
of er and p the simplifying assumptions and the results are consistent. 

In equation (17), then, we replace all square roots by | (3 | and 2 — O 2 
by s , obtaining 



:.2 a 2 , 2 , A + tanh | | x 1 — B 



with 



and 



- 0o cot' + 



A = 



A + 1 " 1 - B tanh | | Bo ' 

Mcff/Mo M + « Sgll 



1 — p H sgn Mo 



5 = -e/eo 



where sgn = 1 for > and sgn = — 1 for < 0. Taking first the 
simplest case of no separation (x = 0), this becomes 

, = 0p 2 (l + e/eo) cot 2 g = 0p 2 (l + e/ep) cot 2 j > 
1 - p g sgn j, mo 



Meff/Mo M + K Sgll 

Substituting the expressions (6) and (1) for p e ff/Mo and p ff we are led to* 



+ - 0o cot *• |/1 (1 + e/eo) )/ ^TT$ 



(18a) 



0_ = -0o cot *• j/I (1 + e/eo) |/ ff °"_ l+ + p/2 (18b) 

* Since reversal of the magnetization is equivalent to interchange of the 
propagation directions, we are at liberty to consider a and p always positive, and 
to deal with the two cases /3 >0 and < separately. 



■ 
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The 0+ mode propagates for all <r, °" . — iz declining steadily from 

to unity. The /3_ mode on the other hand is cut off, with /3 2 = 



1 4- p/2 

at a = 1 — p and then reappears with /3" = qo at a = 1 — p/2. The 
behavior of the two modes is shown in Fig. 5. Self-consistency requires 
that /3 2 ~» co 2 e/io or w 2 e | /Xeff |, whichever is the greater. The condition 
/3 2 » co 2 « I Meff | fails to be met near 



- « - |/f 



*-i-f. 



when 



' T " "~ < 



p Vp-- ± /p _ \ (1 + e/eo) cot 4> 

The condition /3 2 » « 2 eju will fail for /3_ near a = 1 - p. The range for 
this to occur is given by 

P/2 ' 



1 - p - o- < 



1 (1 + e/e ) cot * - 1 



The extension of the Polder formulae to the case of a lossy ferrite 
was given in Part I, (Section 2.1). From the results given there one 
may write 

, a , . _ o"(l + « 2 ) — sgn + ja 

H + k sgn = 1 + v 2n i 2 \ 1 ■ 9o - 

where a is a damping parameter. Substitution of these expressions in 
the slow wave formula for /3 2 will give the effect of loss on the propagation 
constant. In Fig. 6 the complex value of /3_ is shown for a = 0.1 and 
several values of p. The imaginary part of /3+ is small and varies only 
slightly with <r. Fig. 7 shows the initial loss (<r = 0) for three values of a 
and a range of p values. From a knowledge of as a function of a = q} h /u 
and p = o)j//o) it is possible to calculate the loss, Im /3_/(/3 Vl + «/«o 
cot ^), as a function of frequency when the magnetic field and satura- 
tion magnetization of the ferrite are held constant. Fig. 8 shows the re- 
sults of such calculations. It should be noted that the horizontal scale is 
linear in a or 1/w and that the vertical scale implicitly contains the fre- 
quency in the form of l//3o . Both of these distortions of scale tend to 
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Fig. 5 — The non-reciprocal propagation constants of the flat helix. The dotted 
curves represent /3 + , the solid curves /3_ . (Loss free case). 

produce an appearance of sharpness in the variation of the loss at 
higher frequencies. 

If the .slow wave assumption be made again it is possible to obtain 
solutions for the case in which the ferrite does not touch the helix and 
the latter is lossless. With the slow wave approximation, (17) becomes 



& 1 



1 + 



M + k sgn 



to 



Mo 



+ tanh I /3 I x 



/So* COt 2 * e , , . 

1 + — tanh 1/3 I Xo 

CO 



H + k sgn /3 



Mo 



+ 1 



If we write | | .r = v, then the above equation expresses Q± in terms 
of ii. At the same time x = u/\ 8±(ii) | and we evidently have a para- 
metric representation of /3 ± and .r . The results of such computations 
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Fie. 6 — Real and imaginary parts of the reverse propagation constant /3_ 
of a flat helix versus a for various p and for a = 0.1 (the parameter a is marked 
along the curves). The forward propagation constant has a very small imaginary 
part which hardly varies with a. 

are shown in Figs. 9(a), (b) and (c), where is plotted against x for 
various fixed a for two values of -p. It is to be noted that the introduction 
of any characteristic length or scale into the problem, such as is provided 
here by the distance xa immediately produces a great complication in the 
mode spectrum. The plane helix with the ferrite in contact may be 
thought of as a highly degenerate problem. 

To carry out loss calculations using the appropriate expressions for 
M + k sgn /3 would be very tedious in the separated case. However, it 
was pointed out in Part 1 that to order a the expressions for y. and k 
are given correctly if we put a + ja in place of a in the lossless formulae 
and that, in consequence, for small a, the imaginary part of the propa- 
gation constant is approximately given by 
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Figs. 10(a) and 10(b) the loss calculated in this way for the cases 
isidered in Figs. 9 is shown. 
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3. CYLINDRICAL GEOMETRY 

3.1. Impedances 

In this section we consider systems with cylindrical symmetry about 
the propagation direction. All boundaries, those of the circuit as those 
of the medium, are coaxial circular cylinders, and the medium is as- 
sumed to be magnetized circumferentially. The practical means for 
bringing about such a magnetization — for example, thin wires threaded 
through a cylinder of ferrite and carrying a dc current, Fig. 1 (b) — are 
assumed to effect the electromagnetic field to a negligible extent. As 
in the case of planar geometry, we restrict ourselves to fields which have 
no variation along the magnetizing field; that is, in the azimuthal direc- 
tion. Only the ferrite is considered here; the results for a plasma are 
obvious corollaries. The magnetizing field and the dc magnetization 
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Fig. 7 — Real and imaginary parts of & + and /3_ for a flat helix at the very small 
magnetizing field required to saturate the ferrite. 
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Fig. 8 — Attenuation of reverse wave versus frequency at a fixed magnetic 
field and saturation magnetization for various values of w/j/wm = {y.oH)/M. 
The curves for 0.5 and 1 are discontinued when p reaches unity. 

are supposed to be independent of radial distance from the cylinder axis. 
For the geometries employed in practice, this will be a reasonably good 
assumption. Thus it is possible to relate the components of B and H 
in cylindrical coordinates (r, <p, z) through the tensor 



M 





-JK 





Mo 





JK 





M 



where /z and k are given by the Polder relations (1) and are independent 
of position. 

Written in cylindrical coordinates, Maxwell's equations in the ferrite 
are therefore 
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(20a) 
(20b) 



j0H, = jweiE r , 
dH z 



or 

1 a 

- — rH r , = jo)eiE z , 
r or 

tfE 9 = -jMH r - j PH H z ), 

-jfiEr- ^ = ->/*#„, 
Of 

1 !\ 

- — rE v = -juniJPnHr + H,). 
r dr 



(20c) 
(20d) 
(20e) 

(20f) 




1.3 1.4 1.5 



Fig. 9 — Propagation constants for a flat helix separated from the ferrite by a 
distance xo for various values of a. (Loss-free case) (a) /3_ and/3+ for p = 0.2, (b) 
P- for p = 0.8, (c) p + for p = 0.8. In (a), above, the dotted lines bound all <r values. 
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Fig. 9(b) — See Fig. 9. 



As in the case of planar geometry, the field-components can be grouped 
into TE and TM parts; only the TE part will depend explicitly on the 
gyromagnetic properties of the medium. Equations (20a, 20c) and (20e) 
determine the TM field. E T and H v can be eliminated from them, yield- 
ing the familiar wave equation 



r dr dr 



(21) 
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Fig. 9(c) — See Fig. 9. 

whose solutions are zero order Bessel functions (or linear combinations 
thereof) of a kind depending on the region under consideration. Thus 
if j8 > ft , and the region occupied by the medium includes the cylinder 
axis, the modified Bessel function I (ry/(3 2 — ft 2 ) has to be chosen if 
the field is to be finite at r = 0; if the medium extends from a finite r 
to infinity, the function K (r\//3 2 — ft 2 ), regular at infinity, is se- 
lected. Corresp ondingly, if /3 < ft , the Bessel and Hankel function 
-A) (r-v/J8i 2 - /3 2 ) and Ho 1, ' (r V^i 2 - 2 ) replace I and K respectively. 
In terms of the appropriate solution of (21), the remaining field com- 
ponents are 



/Si 2 - 2 dr 



H v = — 



jcoei dEz 
ft 2 - /3 2 ~dr 



The tangential admittance for the TM field is thus 



v _ f — 



jug 

ft 2 - /3 2 dr 



\ozE z 



(22) 
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Fig. 10 — Ratios of reverse attenuation for a flat helix separated from the 
ferrite by a distance x to the reverse attenuation at x = 0. Computations made 
from the approximate slope-formula for the loss, where applicable, (a) p = 0.2. 
(b) p = 0.8. In Ca) all applicable a values lie in the shaded region. 

For example, if /3 > ft , and the medium occupies all space from a finite 
r to infinity 

,7'a>ei K o(c*ir) _ —jcceiKi(air) 



*TM — 



(23) 



ai K Q {air) a\Ka{a\r) 

where a y = V/3 2 - ft 2 . 

The field components of the TE field are determined by equations 
(20b), (20d) and (20f). Elimination of E v and H T from these gives 

!* r ^' + (V- / j'_eZ0V = O, (24) 

r dr or \ r / 

where /?/ = co 2 /x«i (1 — Pa 2 )- The term /3/ 2 — /3 2 in the bracket is already 
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familiar from the planar case; it depends on the magnetic field and on 
the propagation direction in a purely reciprocal way. The term pn(3/r, 
however, reverses sign when either magnetizing field or propagation 
constant changes sign. The solutions of equation (24) are therefore 
different for opposite propagation directions (or opposite magnetiza- 
tions). Thus, in contrast to the planar case, where non-reciprocity arose 
only through the boundary conditions, the cylindrical case is inherently 
non-reciprocal. 

In the absence of the last term in the bracket, equation (24) would 
be solved by zero order Bessel functions, just like equation (21). In the 
presence of this term, the solutions become confluent hypergeometric 
functions. Different changes of variable bring these functions into forms 
known by different names and notations. One such change leads to 
Laguerre functions, another to Whittaker functions. We shall choose 
the latter representation, since it is closely related to Bessel functions, 
and numerical tables seem about equally scarce for all representations. 
In equation (24), let /3 2 — /3/ 2 = a 2 2 , and let a 2 r = y. Then it becomes 

I'.^-Yl-feV-ft (25) 

y dy dy \ y ) 

where x = —@pn/2a2 . Further, let y = x/2, and H t (y) = h(x)/\/x; 
then equation (25) becomes 

d 2 h . / 1 , x l v 



which is the standard form of the equation for zero order Whittaker 
functions. It is a special case of the equation for /z l order Whittaker 
functions : 

^+(i^+^_iK = (27) 

dx> + \ z 2 + x 4/ /l °' 

The solution of equation (27) which is regular near zero is denoted in 
the literature by M x ,„(x); it is proportional, in the limit x = 0, to the 
Bessel function 7 M (x/2). The solution of equation (27) regular at infinity 
is denoted by W Xtli (x), and in the limit x = 0, is proportional to K p (x/2). 
The factors of proportionality are found in Appendix I. 
In this notation, the solutions for H z are thus 

M x . (2a 2 r) W Xt0 (2a 2 r) 4 



* If P < Pt , both argument and suffix x are imaginary. These functions are 
then related to J and Ho respectively. 
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Once the appropriate H z is determined, H r and E v are given by 

Br - m ' (f> 9 -£- »'l***B.) , 

/3 2 — w 2 juci \ dr . / 

/3 2 — «Vi \ dr 
The tangential impedance for TE-fields is thus 



(28) 



Ef, w/i 



(iloBff.-prf). 



The reader will recall that for isotropic media (p fl = 0) the numerator 
of the right hand side of equation (29) can always be expressed as the 
ratio of first order to zero order Bessel functions by virtue of relations 
like Kq (x) = -Ki(x), lobe) = h(x) and so on. Analogous results 
hold true in the present case. Suppose, for example, that Ave are dealing 
with a geometry such that the correct H z is 

W x , Q (2a 2 r) 
H z = — yjr — = R x0 (2a 2 r), say 



Then 



and 



— — — = 2ot2R x o/R x o , 
rii or 



2a 2 "/i i2 x o + X#*o 



i(/3 2 - a>Vi) # x o 

It is shown in Appendix I that 

Rxo' + X#xo = (X - I) fi x i • 

Therefore 

_ co>ja 2 (2x — 1) R x i(2a 2 r) 
TE iCS 2 - coVeO R x0 (2a 2 r) ' 

= <W2 X - 1) TF Xil (2a 2 r) 
J'03 2 - coVO W xfi {2a 2 r) ' 

A similar difference relation shows that if the region is such that 

H t = M x , (2a 2 r)/\/W, 



(30) 



(33) 
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then 

7 = fr>M«2(3^ - x)M Xil (2a 2 r) , ^ 

TE j(f3> - coVe 1 )M x .o(2a 2 r) ' K * 1} 

In the unmagnetized case equations (30) and (31) reduce to 

_ o>/zoq:i Kijair) _ _ mxoif i(air) , . 

m - coVoei) K ( ai r) jotiKoM ' ^ ' 

and to 

w/zo h(ctir) 
jai 7 (o!ir) " 

One might be led to believe that the search for solutions of Maxwell's 
equations with angular dependence e 3 "* will lead to n th order Whittaker 
functions (just as in the isotropic case this dependence leads to n th 
order Bessel functions). Such is not the case, however. Unless n = 0, 
one is led to two simultaneous second order equations for E, and H z , 
and the character of the problem is changed completely. 

3.2. The Cylindrical Helix 

We are now in a position to derive the characteristic equation for a 
close wound cylindrical helix and approximated by a helically conducting 
sheet surrounded by ferrite. We confine the discussion to the case in 
which the ferrite is in actual contact with the helix, Fig. 11; the case of 
finite separation discussed for the planar helix (Section 2.3) would be 
too cumbersome here. Losses are neglected. If they are not excessive, they 
can be deduced from the curves for the propagation constant in the loss 
free sample by differentiation, as outlined in Sections 2.1 and 4.15, 
Part I. 

The boundary conditions are just the same as in the planar case. In 
Section 2.3 they were stated in terms of admittances, and it is only 
necessary to substitute for these the admittances just derived for cy- 
lindrical geometry. Thus for H y /E z we substitute Ftm , and for H z /E v 
we write F T e = 1/^te . 

If superfixes i and e refer to the inside and outside of the helix (in 
Section 2.3 on the plane helix i and e were denoted by " — ", and "+"), 
the characteristic equation is 

{Y TM li) - F TM (e) U cotV = (r TE (0 - F TE (e) ) r=ro , (34) 

where r is the radius of the helix. 



* The appearance of different factors (2\ — 1) and (34 — x 2 ) is simply due to 
the way the functions W, M are normalized in the literature. 
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For waves bound to the helix, Yte M is to be derived from that solu- 
tion of (24) which tends to zero ar r — > =©. This solution is W xfi (2a 2 r)/ 
■\/2a 2 r, and so FxE (e) is given by equation (30) 






1 = j /3 2 - aw W x , (2a 2 r) 
Z TE com a 2 (2 x - 1) W x ,i(2a 2 r) ' 



Similarly, from equation (23), we have, for bound waves, with E z ~ K Q , 



FERRITE 




DIRECTION OF 
CONDUCTION 




CIRCUMFERENTIALLY 
MAGNETIZED 



Fig. 11 — (a) Cylindrical helix surrounded by ferrite, (b) Magnetic field lines 
projected onto a plane containing the axis of the helix. 
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(«) juei Ki(air) 



Y. 



TM 



a x Ko(ctir) ' 



Inside the helix, we require the solution for an isotropic region which 
remains regular as r — > 0. Accordingly 

v (») _ i«o h(a r) r- — 

*TE - Y~< \ > a ° V P* - ft> 2 > 

and 



Q 2 2 

Pq = CO fl € 



v (,) _ • w «o /i(«o?"o) 
a lo{a ro) 

where eo , Mo are the dielectric constant, and permeability of vacuum. 
Combining these expressions in (34), we obtain after slight rearrange- 
ment 



h(aoro) , «o ci Ki(ociro) 



2 I (oioro) Moao/3 — w M«i W x , (2a 2 r) 
«o 



(35) 
tan 2 ^ 



Io(oc r ) ai€ Ko(ociro) _ 7i(a r ) m«2 (2x — 1) W x ,i(2a 2 r)_ 

which determines /3. 

A complete solution of equation (35) is out of the question. However, 
as in the planar case, for the slow waves used in traveling wave tube 
work, the equation may be simplified so that solutions may be com- 
puted rather easily. For electron velocities usually employed the result- 
ant /3 must be about 10/3o . Therefore in equation (35) it will be permis- 
sible to neglect all the quantities /3 ", ft , /3a , co /iei , in comparison with 
2 , except in the narrow ranges of magnetic field such that y. or 
/x(l — pu) becomes very large. This will occur near ±o- where <t = 
— p/2 -f Vp 2 /4 + 1> ano - near o" = 1. A solution obtained by as- 
suming a large /3 must be self-consistent; that is, it can be credited only 
in regions where it does, in fact, predict large 0. However, in Section 
2.3 it was shown for the plane helix that in any practical case the ranges 
of magnetic fields so excluded are very narrow, even in the loss-free 
case, and one may suppose that this is true also in the cylindrical case. 

For slow waves, each of the a's reduces to | /? | ; the absolute value 
sign derives from the fact that the positive square root is implied in the 
definition of the a's. Therefore 

*■+¥$- ''I***?. (36) 

Now the suffix x of the Whittaker functions no longer depends on the 






972 



THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1954 



magnitude of j3, and it is chiefly for this reason that further progress is 
possible. For large 0, equation (35) can be written 

JiQfllro) 6i gi(l/3|r ) 
Jo(| | r ) 60 g (| g | r„) 

1 PF x .o(2 | ^ |r ) " (37) 



/3 2 = ft 2 cot 2 ^ 



/o(l g 1 r ) _ 

/i(| /3 | r ) m. (2x _ 1} TF x .o(2 | p | r ) 

Mo 

where x is now given by equation (36). Equation (37) is now solved by 

-7 




0.2 0.4 



1.2 



0.6 0.8 I.* 

r /3 COT V 

Fig. 12 — Reduced forward and reverse propagation constants versus reduced 
radius of a cylindrical helix (loss-free case) for various a and p. The range o\ 
< a < <r where 

contains an infinity of shape resonances and is not shown here, (a), above, p = 0.2. 
(b) p = 0.6. (c) p = 1.0. 
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Fig. 12(b) — See Fig. 12. 



1.2 



1.4 



1.6 



the following procedure: Introduce the parameter 



= I I n 



(38) 



For a given p H and u (or <r and p), and a given sign of /3, each value of 
u determines /3 through equation (37), and then r through equation (38). 
Thus /3 can be plotted versus r . The procedure is repeated for the op- 
posite sign of /3 (and therefore the opposite sign of x). A different curve 
of jS versus r is then obtained. Thus for a given value of r , the "forward" 
and "backward" propagation constants are different in magnitude. 
The results (computed for a typical ratio ei/« = 10) are conveniently 
stated in terms of /3 = 0/(8o cot \f/) and f = r /3 cot \J/ and are shown 
in Fig. 12(a) to (c), and again, for fixed f , in Fig. 13(a) to (e). We note 
that for fo in excess of about 1.5, the results are almost the same as those 
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Fig. 12(c) — See Fig. 12. 
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Fig. 13 — Reduced forward and reverse propagation constants versus a for 
various reduced radii of a cylindrical helix. (Loss free case). The region 
1 - p < a < ff is omitted, (a), above, p = 0.2. (b) p = 0.4. (c) p = 0.6. (d) p 
= 0.8. (e) p = 1.0. 
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for a flat helix. This is due to the fact that the large dielectric constant 
reduces the circuit wavelength so much that the radius appears infinite 
by comparison once it exceeds 1.5. In traveling-wave tube practice, 
however, f is generally below 1.5. 

The behavior of the /3 — f a curves can be understood from the be- 
havior of the ratio 

Z X &1 = zm 
W X .M 

and of the coefficient 1/[m/mo (2x - 1)] • Suppose first that x is positive. 
When x exceeds zero only slightly, Z x (u) behaves essentially like 
Kq(u/2)/Ki(u/2). This function is always positive; it begins at when 
u = with a vertical tangent and steadily increases to unity asw-» oo . 
Z x (u) varies in the same way, see Figs. 14 and 15, in the range < x < 
Y%. For Y 2 > x ^ /-i,W x ,o, and therefore Z x , has a zero which increases 
from u = at x = V A to u = 1 at x ■ %• Accordingly Z x (u) in % > 



2 x 2 




Fig. 14 — Schematic behavior of the function Z x (u) = W x , a (u)/W x . l (u) in 
various ranges of x- 
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Fig. 15 — The function Z x (u) versus x for various u, in the range — % <x 
< H- 

X> }4 starts from at u = with a downward-directed vertical tangent, 
achieves a negative minimum, then increases, through its zero, to the 
asymptotic value unity asM-> oo. The minimum becomes deeper as 
X approaches %. At % = % % W Xil (u) developes a zero at u = 0, which 
steadily moves to larger u as x increases further towards %. At the 
same time the zero of W Xi0 already discussed moves from 1 to 2 + \/2, 
and a new zero arises at u = 0, x = %, which increases to 2 — y/2 
as x approaches %, but which lags behind the zero of W x ,i{u). The 
function Z x (u) now has a pole and two zeros, and behaves as shown in 
Fig. 14. This process continues; each time x passes (2n + l)/2, a 
new zero and a new pole appear. (For a detailed list of poles and zeros 
the reader is referred to Appendix I). To apply these results, we first 
resort to the Polder relations. 

In terms of a, p, we have, for /3 negative 



and 



x = y 2 



V 



1 — pa — 



1 - 



*-(2 x -l) P-! + * 
Mo 



= A, say. 
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The characteristic equation is 
7i(w/2) 
/3 2 = 



€iKi(u/2) 



7o(t*/2) 6qKq(u/2) . |^i r - = w / 2 



(37) 



and can now be discussed in terms of a at a fixed p. Suppose that p < 1. 
Then 4 is negative for er < 1 — p. In the same range, x < H, s0 tnat 
Z behaves essentially as K /Ki . Therefore both numerator and denomi- 
nator are positive for all u, and the ratio tends to the planar result 

l+a 

p ~t=a 

as u _» oo . As u — > 0, j§ 2 tends to zero along the vertical, as can be shown 
by an examination of the various functions near u = 0. For a < 1 — p, 
the course of the j8 2 versus w-curves is as shown schematically in Fig. 
16(a), and it is easily seen that the |8 versus f curves run in essentially 
the same way, Fig. 12(a) to (c). However, as a approaches 1 - p, the 
iS 2 versus u curves steadily fall, until at a = 1 — p, J3 2 = for all finite 
u, since A ■= — » . 

* As o- passes 1 — p, A changes sign and at the same time x passes ^ 
so that Z x acquires a zero. As a varies from 1 - p to 1 — p/2, A decreases 
from + oo to unity. Therefore, while u < u x , the zero of Z x , 1 — AZ X 
is positive; however as u increases beyond wi , {h(u/2)/Ii{u/2)) — 
AZ x (u) eventually passes zero, since Z x (u) and T Q (u/2)/Ii(u/2) both 
approach unity. On the other hand the numerator of equation (37) is 
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Fig. 16 — Schematic variation of /3_ with u. a) x < ^5 b) ^ < x < M; c) % 
< x < %. 
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always positive; therefore /3 s approaches infinity at 7 (w/2)//i(w/2) — 
AZ x (u) = 0, and no real values of /3 exist thereafter (see Fig. 16b). 
Since this "cut-off" occurs at a finite value of u, the corresponding value 
of r is zero. This explains the bulging of the corresponding j8 — fo curves 
in Fig. 12(a) to 12(c). 

The next major change in the curves occurs when x exceeds %, (that 
is, a exceeds 



tri 



2 T 



v- 



i + ?-S 



For p < 2, o-i is still less than 1 — (p/2), so that, initially at any rate, 
A is still greater than unity. In addition to the infinity of /3 2 just dis- 
cussed, a further infinity arises between u = 0, and the pole of Z x (u), 
as is seen from Fig. 16(c). /3 2 increases from zero at u = to this infinity, 
thereafter it is negative, until the pole of Z x (u) is reached. Thereupon 
it resumes at /3 = and approaches infinity at the zero of the denomina- 
tor [Io(u/2)/Ii(u/2)] — AZ x (u) already discussed. Thus there are now 
two branches of the /3 2 — u curve; their corresponding traces in the 
/3 — f plane are shown schematically in Fig. 17. (The computations on 
which Fig. 12(a) to 12(c) were based were broken off at a = a x .) 

A further branch is added each time x increases beyond a number of 
the form (2n + l)/2 (a increases beyond 



= -£ + 

2 T 



/• 



+?- 



V 



2n + 1 



These all resemble the two branches just discussed, until n > 




Fig. 17 — Schematic variation of /3_ with r for % < x < H- 
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Fig. 18 — (a) /3_ versus u when 1 — (p/2) < a < <r . (b) /3_ versus r under the 
same circumstances. 



Y % (- - 1). When this occurs, <r„ > 1 - (p/2), so that there will be a 

P 
value o- between a n -\ and <r„ beyond which the denominator no longer 

decreases through zero as u — > w , but approaches a finite positive value, 
Fig. 18(a). Accordingly /3 2 approaches a finite positive value, and cut-off 
of the extreme right hand branch, Fig. 16(c), no longer occurs. The cor- 
responding |Q versus fo branch is as in Fig. 18(b). 

As n — > oo (o- n — > ao = Vl + (p 2 / 4 ) — p/2) the number of branches 
increases to infinity. This situation resembles that in the completely 
filled waveguide (Part I), where we found an infinity of modes ("Shape- 
resonances") in the range <r < <r < 1. In the present case, however, 
they are to be found in the range 1 — p < c < <x . 

When o- = o-o + 0, x is infinite and negative. The function Z x (u) is 
then constant and equal to unity. A is less than unity, and the denomi- 
nator of equation (27) has no zeros. The /3 versus f curve is now "nor- 
mal" again, see Fig. 12(a). As a increases further, the curve falls (since 
A decreases steadily to — 1 as a — * °o), and no more qualitative changes 
occur. 

3.3 Cylindrical Waveguides 

As pointed out before, the fact that the propagation problem in the 
cylindrical case can always be integrated in terms of Whittaker functions 
when the fields show no angular variation is an accident, and in view of 
the lack of numerical tables, not a particularly fortunate one. Only in 
special cases (like that of the slow- wave helix) is the text-book informa- 
tion on these functions of any great utility. In general, it will be more 
convenient to solve the differential equations numerically. However, 
for completeness, we shall state some of the formal results for a cylin- 
drical waveguide containing a cylinder of circumferentially magnetized 
ferrite, and propagating a TE mode. 
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First we consider a waveguide, radius r , into which is fitted a hollow 
cylinder of ferrite, outer radius r , inner radius n . In that cylinder, the 
magnetic field H z may be taken to be a superposition AS x0 (2a 2 r) + 
BR x0 (2a 2 r) where, as before, 

a \ = /3 2 - co 2 e lM (l - p„ 2 ); S(z) = —7^ ; R(x) = 



Vi ' rtW " Vi 



2a 2 



a 2 may be either real or positive imaginary. [In choosing this combination 
we depart from the usual practice of taking a superposition of J and 
No in the isotropic case. Were we to follow this practice, it would be 
necessary to define a new function R xll (2jx)e 3(,l " l2) + R xil ( — 2jx)e~"' r 
to correspond to N?(x). Our choice corresponds to taking a combination 
of Jo(x) and one of the Hankel functions H (x) in the isotropic case. 
Since the functions H, J, N are linearly dependent, this will not affect 
the results.] 

In view of the difference relations, equation (39) in Appendix I, 
and of equation (29) we obtain for the impedance in the ferrite 

§, m jWMgg [A(K - x 2 )S xl (2a 2 r) + B(2 X - l)R xl (2a 2 r)] 

H t OS 2 - o>V0 [AS x0 (2a 2 r) + BR x0 (2a 2 r)] 

A and B must be adjusted so that this quantity vanishes at r , the guide 
wall. This gives 

$ = -<2x - Wi - x 2 ) **""• 



H x /3 2 — co 2 /x«i 

W x .i(2a2ro)M x ,i(2a 2 r) - M x , 1 (2a 2 r )ir x , 1 (2Q !2 r) 
(2 x _ 1 )Tr Xll (2a S r )M x .o(2a^) - Q£ - X 2 )M x , 1 (2a x r )W x , (2a 2 r) ' 

In the vacuum, between r = and r = ri , H, is Io(a r) and the im- 
pedance is EJH, = —j(wno/ao)Ii(a r)/Io(a Q r) where a 2 = /3 2 — w 2 m«o . 
At r = r x , the ferrite-vacuum interface, the two impedances must be 
equal. Thus we obtain the characteristic equation 

- rrvi " (2 * - 1)( ^ - **> stt- 

IF x , 1 (2a 2 ro)M' x , 1 (2a 2 r 1 ) - ■M" x .i(2a 2 r n )IF x , 1 (2 a2 r 1 ) 

(2 X - l)T^ Xil (2a 2 r )ilf x ,o(2a 2 r 1 ) - Q4 - x 2 )M xA {2a 2 r )W x . a (2a 2 n) ' 

(It is understood that for "normal" waveguide propagation a will be 
imaginary, and the I will be replaced by J). As a simple illustration we 
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consider the case in which the ferrite cylinder fills the waveguide com- 
pletely. E 9 is then proportional to M x i(2a 2 r), and so is determined from 
the condition 

M x i(2a 2 r ) = 0. 

When "normal" propagation pre vails (ff less than the natural propaga- 
tion constant of the medium w-y/Wl - P// 2 )) Dotn a * and X are imagi- 
nary, equal to ja% and jx', say. Under these circumstances little is known 
about the zeros of M. However, it is possible to say something about 
the solution for large radial mode numbers. It follows from Erdelyi 
et al. 2 1, p. 278, formula (2), that for large argument 

M ix >,x{2ja 2 'r ) = const • sin [a 2 'r + x'log a 2 'r + x'log 2 + $ (%') - */4]. 

where $(x') = arg r(% + jx')- The zeros of this expression are at 

a 2 'r + x log a 2 'r = ~ - x ' log 2 - *(x'). w = a large integer 

This equation may be solved graphically by setting a 2 'r Q = u, assigning 
values to p, p B , u (and hence to x', ct 2 ). From a solution u one then 
finds 

r = u/a 2 . 

M also has zeros for real a 2 , x, if X is large enough. Thus the wave- 
guide will support waves with a Z3 2 greater than oj/*€i(1 — p H ). It is 
shown in Reference 2 (1, p. 289) that when x is between K and ^, M 
has one zero, when x is between % and %, M has two zeros and so on. 
Suppose that p n is negative, = — \ p H \ • Then 

= gj_ggj 
X V/3 2 - /3 2 2 " 

For real positive /3, this equation has a solution for if | pa \ < X- 



J8 = 



Vx 2 - PB 



If M < X < Yi, Af will have a zero w(x) depending on the value of x- 
Thus the equations 

„ X& "OP _ VX 2 - PR 2 ,.(„\ 

t> = vx-p/ ro = vf=j* \ PB \fr u{x) 

solve the propagation problem parametrically. Similarly when x is 
between % and %, there are two zeros of M given by two functions 
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Mi(x) and W2(x)- There are now two possible modes, with the same 
restrictions on p„ . An additional mode arises each time x is allowed to 
pass a number of the form (2n + l)/2. It is to be noted that these 
modes are not confined to the resonance range. For /3 positive, they can 
exist in the range co > <r > a and in the range — o- < a < 0. 

Appendix I. Some Properties of Whittaker Functions Used in 

This Paper 



I. RELATION TO BESSEL FUNCTIONS 

Lt M^m _ 2 ,„ r0j + 



u W x „(2x) K,{x) 

x-o 



V2x Vx ' 



, W Xtl l ( 2jx ) -y/^ 

£ V^i * "2" r *' W ' and 

x ^o|_ V2ja; V-2jx J 

II. DIFFERENCE RELATIONS 

The following results can be obtained either by reference to Erd61yi, 2 , 1 
pp. 258, 254, by differentiation and subsequent integration by parts 
of integrals such as 

u+l/2 -x/2 f oo 

W X M = X * , e^^Wl + r)" + - 1/2 dr 

T(m + Yi - x) Jo 
or by observing that combinations of the form 
/- d W x ,o(x) 

^ x Tx Vx +*"*#> 

satisfy Whittaker 's equation with y. — 1. In the last mentioned method, 
the required constant multiplying the first order Whittaker function 
can be obtained by reference to the limiting behavior for small x. If 
R X n = W x ,p(x) / -\/ x and S xll = M Xili (x)/\^x the results are 

RW + xR x o = (x - V 2 )R x i , , v 

(39) 
>V + xS x0 =H(M~ x) S xl . 
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III. ZEROS 

When x = (2n + 1)A n = 1, 2 , W x .„/x?~ 112 reduces to a 

polynomial times a function of x- This may be inferred from the asymp- 
totic expansion, 



(l+ t 

\ n-l 



lV - (x - y 2 )V - (x - ^) 2 ]---[M 2 -(x-n + ^) 2 ] 



n!z n 



)■ 



which terminates if n = 1 and x = n + }4 ( n = 1, 2 . . . ), or from the 
fact that for these values of the suffixes, W reduces to the generalized 
Laguerre polynomial 

Similarly when x = (2n + 1)A n = 0, 1, 2 , TT X|0 reduces to the 

Laguerre polynomial 

L x - {m (x). 

The zeros at the critical values of x are given in the following table 



X 


Zeros of I^x.o 


Zeros of W x , i 


V, 









None 


■¥, 


0; 1 









V, 


0; 0.586; 3.414 






0;3 


% 


0; 0.416; 2.294; 6.290 






0;6;2 


% 


0; 0.323; 1.746; 4.537; 9.395 




0; 1.517; 4.312; 9.171 


m 


0; 0.26356; 1.413; 3.596 


7.086; 


12.641 


0; 1.227; 3.413; 6.903; 12.458 



Between n + K and n + %, W x , has n + 1 zeros (n - 0, 1, 2 . . .) 
and TF x ,i has n zeros (n = 1, 2 . . . .). 

The zeros of ilf x ,i coincide with those of W x ,i when x = n + /4, and 
at those values of x only. The functions M Xli and W xA then are propor- 
tional to each other.* 

IV. THE RICCATI-EQUATION FOR THE IMPEDANCE FUNCTION W X ,JW X ,\. 

The computations concerning the cylindrical helix required a study 
of the function 

7(oA - W xM 



* At these critical values of x, the solution to the problem of the hollow cyl- 
inder of ferrite in the waveguide breaks down, since M and W are then not in- 
dependent. A further independent solution must then be constructed. 
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We will show that Z x (u) satisfies a non-linear first order differential 
equation of Riccati-type. From the difference relations, equation (36), 
we have 

( Y _ w\ w xM m Wj(u) _ 1 

W Xt0 (u) " W x , (u) 2u T * 



say 



and from Whittaker's equation 

d (WJ\ (WJ\ ( ] 
du\W x ,J^\W xfi ) T \4m 



4w 2 w 



-.).« 



Therefore 



or 



Finally, let 






i-i+S-^f+tf-K)^. 



£ = (x - V 2 ) g(u) = 



W X M 
W xA (u) ' 



Then 



£ifc-» + 



(i - 2 x )z + ( x + K) Z 2 . 



Since this equation is satisfied by M Xi0 (u)/M x ,i(u)) as well as by 
W Xi o(u)/W Xt i(u), a selection has to be made from all the possible solu- 
tions of this equation. We require the one which for large u approaches 
unity. But for large u the equation is 



^ = (i-z)[ x -y 2 -( x + y 2 )z], 



whose integral is 



Z = 



(x - V2)Ae u - 1 



(x + MMe" - 1 

For large u, therefore, the solution is either unity, when A = 0, or else 
(x — H)/(x + %), A z}z 0. The solution with A =£ corresponds to the 
M functions; that with A = to the PF-functions. 
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The case A = was integrated on an analogue-computer, and the 
results are shown in Fig. 14(b). The computation was restricted to the 
range | x I < M- Beyond these values, the helix-problem was discussed 
only qualitatively. 
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